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ABSTRACT 

Active polarimetric imagery is a powerful tool for access- 
ing the information present in a scene. Indeed, the polari- 
metric images obtained can reveal polarizing properties of 
the objects that are not avalaible using conventional imag- 
ing systems |1, 2|. However, when coherent light is used 
to illuminate the scene, the images are degraded by speckle 
noise. The polarization properties of a scene are character- 
ized by the degree of polarization. In standard polarimetric 
imagery system, four intensity images are needed to estimate 
this degree |3|. If we assume the uncorrelation of the mea- 
surements, this number can be decreased to two images using 
the Orthogonal State Conti-ast Image (OSCI) |4|. However, 
this approach appears too restrictive in some cases. We thus 
propose in this paper a new statistical parametric method to 
estimate the degree of polarization assuming correlated mea- 
surements with only two intensity images. The estimators 
obtained from four images, from the OSCI and from the pro- 
posed method, are compared using simulated polarimetric 
data degraded by speckle noise. 

1. INTRODUCTION 

Polarimetric imagery consists in forming an image of the 
state of polarization of the light backscattered by a scene. We 
consider in this paper that the scene is artificially illuminated 
with coherent light (laser). For example, this illumination 
is used in active imagery in order to combine night vision 
capability and to improve image resolution for a given aper- 
ture size. In practice, using a coherent illumination produces 
speckle noise that deteriorate the image 1 5 1 . However, the 
backscattered light gives information about the capability of 
the scene to polarize or depolarize the emitted light and thus 
allows one to determine the medium that compose the scene. 
These information can be described by a scalar parameter: 
the degree of polarization of light. This quantity is obtained 
in standard configurations of polarimetric systems using 
four pair of angular rotations of both a compensator and a 
polarizer Four transmittance are thus recorded 1 3 1 that lead 
to the estimation of the degree of polarization. However, this 
system is complex and it is interesting to develop methods 
to estimate the degree of polarization that could reduce the 
number of images to register In |4|, the authors proposed to 
estimate the degree of polarization with only two intensity 
images, however this method relies on the assumption that 
the measurements of the two components are uncorrected 
which can be in some cases a too restrictive hypothesis. This 
paper extends the work of (4] by taking into account the 
correlation of the different components. 



Let us first introduce the context of the study. 



2. BACKGROUND 

The electric field of the light at a point of coordinates r (vec- 
tor of 3 component) in a 3D space and at time t can be writ- 
ten, if we assume the light to propagate in a homogeneous 
and isotropic medium, as 



E(r,f) = [A^(r,f)ex+A^(r,f)ey] .e 



-ilnvt 



(1) 



where v is the central frequency of the field and Gx , By are 
unitary orthogonal vectors (in the following bold letters rep- 
resent vectors). 

The terms A^{r,t) and A^(r,r) are complex and define the 
random vector called Jones vector 



A^(r,f) 
A^r,t) 



(2) 



The state of polarization of light corresponds to the prop- 
erties of E(r,f) at a particular point of space. It can be de- 
scribed by the covariance matrix F 



<A^*A^ 



<ArA^ > 



> <A{*A'^ > 



(3) 



where < . > and .* define respectively the statistical average 
and the complex conjugate. For sake of brevity, the following 
notations for F are introduced 
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Cl 
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Let us note that this matrix can be diagonalized since it is 
hermitic. 

In the case of coherent light, the electric field is represented 
by the complex Jones vector A which follows a Gaussian 
circular law |5| 
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where A^ stands for the adjoint of the vector A. 
The degree of polarization is defined by 



(6) 
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where Hi and /X2 are the eigenvalues of F (/ii > /X2 > 0). 
The degree of polarization is a scalar parameter that charac- 
terizes the state of polarization of the light: if f = 0, the light 
is said to be totally depolarized, and if P = 1, the light is said 
to be totally polarized. In the intermediate cases, the light is 
partially polarized. 

With the notations introduced in one can show that (5) 
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The knowledge of this quantity allows one to study the way 
the illuminated scene polarized or depolarized the emitted 
light. This degree gives information about the nature of the 
medium in the scene. In the standard configuration, four 
measurements are needed to estimate it. In the case of un- 
corrected measurements, two intensity images gives a good 
estimation of this degree using the OSCl |4|. However, in 
some cases the uncorrelation of the measurements may be 
not valid. 

In this paper, an original estimation method that both uses a 
pair of images and that accounts for correlation in the compo- 
nents is proposed. We recall in the following the method pro- 
posed in 1 4 1 and we extend it to correlated measurements. We 
then compare them through statistical measures using simu- 
lated data. The results are also presented when the standard 
estimation of is used (four images) as this case is expected 
to give the best results. Finally we conclude and give the per- 
spectives of this work. 

3. THEOSCI 

In the case of two uncorrected images, the Orthogonal State 
Contrast Image (OSCI) is determined from L4J. L6J. L7J 



P(0 



(9) 



where Ii (/) and 12(1) are intensity measurements at the pixel 
site (/), assuming a lexicographic order for the pixels. 
These two images are obtained with simple polarimetric sys- 
tems. First, the scene is illuminated by coherent light with 
single elliptical polarization state. Then the backscattered 
light is analysed in the polarization state parallel (which 
leads to = |Ax(;)P) and orthogonal (which leads to 
^2(0 = I^f(OP) '^o '^hs incident one. 

The OSCI is an estimation of in each pixel provided that 
the materials of the scene modify the degree of polarization 
of incident light without modifying its principal polarized 
state ('). Let us recall that this kind of material is called pure 
depolarizer For such materials, the covariance matrix F is 
diagonal and, since the diagonal terms represent the intensity 
images, the OSCI gives an estimation of P^ with 



A (0-/2(0 
h{i)+h{i) 



(10) 



In the case of non pure depolarizer objects (i.e. F is non 
diagonal), the OSCI still reveals interesting contrast image 

' State represented by the eigenvector associated to eigenvalue fli . 



but no more defines P^. This leads us to a new method that 
considers the cases of correlated measurements. This is the 
object of the following part. 

4. CORRELATED MEASUREMENTS 

In the case of correlated measurements, the covariance ma- 
trix is non diagonal and of the form 0. In its standard esti- 
mation, the degree of polarization needs four measurements, 
however, two images are sufficient to get an estimation of P^ . 
Indeed, the coefficient a 1 is obtained from one measurement, 
as the coefficient 04 and the squared modulus of 02 can be 
estimated from the cross-correlation coefficient 812 between 
two measurements /j and 12- We have 



5l9 = 



hhp{h,h)dhdl2 =< hh > 



(11) 



5i2 can be calculated by using the joined density probability 
function p{Ii,l2) assuming that A is Gaussian circular (i.e. 
the speckle is supposed to be fully developped). It can be 
shown that the correlation coefficient is obtained with 



(12) 



where \c2\ stands for the modulus of C2- Calculating the cen- 
tered correlation coefficient defined by 




A12 = 5i2- </i Xh > 
After some simple algebra, we get 

A12 = |fl2p 



(13) 



(14) 



Thus the coefficient |fl2p can be obtained from two 
measurements with < /1/2 > — < Ii >< I2 >■ This remark 
leads to write the following property. 

Property A: For fully developped speckle fluctuations, 
the degree of polarization can be obtained from only two 
intensity images. 

One can easily note that the degree of polarization can 
be written as a function of the OSCI with 



A12 



(</l > + </2>) 



(15) 



From (IIOL the degree of polarization estimated from the 
OSCI is clearly under-estimated. Thus, we can correct the 
OSCI in order to get an estimation of the degree of polariza- 
tion using the correlation coefficient A12. 
In the following part, the different estimation of the degree of 
polarization are compared to the estimation with four images 
through simulated data. 

5. COMPARISON WITH NUMERICAL 
EXPERIMENTS 

We generated R experiments of samples of complex Jones 
vectors which follow a Gaussian circular law. The covari- 
ance matrix is known and thus P^ is also known. Under the 



assumption that the statistical average can be estimated by 
spatial averages in homogenous regions, the coefficients ai 
can be estimated from a single image ( /i = \Ax\^) like the 
coefficient 04 (I2 = |Ay p) since 



ai 



04 



(16) 



(17) 



where {Ax{i),AY(i)} represents the component of the Jones 
vector for the sample /. 

The estimation of differs in the studied cases by the way 
|fl2p is estimated. Three different methods are used: 
• Case of four images 

In this situation, we have both the real and the imaginary 
part of the coefficient 02- The quantity |fl2p is estimated 
by 



Mi) = 



if A;,(/)A*(0 



1=1 



(18) 



• Case of two images with the OSCI 

In this case a2 is assumed to be equal to zero. 

• Case of two images with the proposed approach 
The coefficient |fl2p is estimated using ( I13t with 

P/(o = ^Lf=,iAx(onAr(op 

(19) 

In the three cases, was estimated from the relation (|8} 
with the estimated parameters. 

In order to characterize the precision of the estimation, one 
considers six examples of matrix F. 
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0.2 + 0.5/ 
6 



;r2^ 
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r3 = 



82 13/ 
-13/ 17 
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18 7 + 8/ 
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5 — 



30 16 - 8/ 
16 + 8/ 14 



6 — 



1.25 
-5.5/ 



5.5/ 
26 



(20) 



(21) 



; (22) 



These matrices were chosen such as the degree of polariza- 
tion are approximatively in {0.2,0.4,0.5,0.6,0.8, 1}. 

The simulations are performed for R = 1000 realisations 
of N = 10000 samples for the six covariance matrices. For 
the two matrices Fi and F5, supplementary cases have been 
studied whenA^ G {100,500, 1000,5000, 10000}. The results 
are presented in figures [0 12] E |3 |5] |6] Several points are 
important to notice. First of all, as expected, the best esti- 
mations of the degree of polarization, regarding all the cases 
tested, was achieved with four images. However the pro- 
posed approach that relies on two correlated images produces 





i 




1 




1 




\ 

1 

\ 



- OSCI 



Covariance matrix 



Figure 1 : The degree of polarization is plotted as a function 
of the six covariance matrix for R = 1000 and — 10000. P^ 
is the true degree of polarization, - A is the degree esti- 
mated from four measurements, f ^ - 1 is the degree obtained 
with the numerical simulations using the proposed method 
for evaluating \a^, P^ - OSCI is the degree estimed from 
the OSCI. 
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Figure 2: The experimental variance of the degree of polar- 
ization is plotted as a function of the six covariance matrix for 
R = 1000 and = 10000. cr - A is the standard deviation of 
the degree of polarization estimated from four measures, a - 
I is the standard deviation of P^ obtained with the numerical 
simulations using the proposed method for evaluating |fl2p, 
(T - OSCI is the standard deviation of P^ estimed from the 
OSCI. 



good estimations of the degree of polarization whatever the 
covariance matrix is used (figCJ as soon as N > 1000 (fig|3l 
and fig|5])- Note that the estimation with the OSCI gives re- 
sults which cannot be used if the term A12/ ((< > + <l2> 
)^) is too high (for example if |fl2p is non negligible). Fig|2 
|4]and|6]show that the experimental variance using four mea- 
surements or the OSCI are comparable whatever F and are, 
whereas the variance obtained with the proposed approach is 
larger However this precision should be sufficient for some 
practical applications. This point should be studied in details 
in a future work. 

6. CONCLUSION 

We have proposed a new approach to estimate the degree 
of polarization on polarimetric images degraded by speckle 
noise. Assuming that the speckle is fully developped, this 
method allows one to estimate this degree with only two in- 
tensity images whereas four images are needed in a standard 
experimental setup. This presents a great interest in term of 
reduction of cost of the imagery system since the original 
setup can be simplified. The proposed approach has been 
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Figure 3: The degree of polarization is plotted as a function 
of the number of sample for R = 1000 for the co variance 
matrix Fi. is the true degree of polarization, - A is 
the degree estimated from four measures, - lis the degree 
obtained with the numerical simulations using the proposed 



method for evaluating |fl2p 
from the OSCI. 
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Figure 4: The product of N and the experimental variance of 
the degree of polarization (a^) is plotted as a function of 
for R — 1000 for the covariance matrix Fi. N<7^ - A is the 
product obtained when P^ is estimated from four measures, 
N<J^ - 1 is the product obtained when P^ is estimated with the 
numerical simulations using the proposed method for evalu- 
ating |a2p, N<J^ - OSCI is the product obtained when P^ is 
estimated from the OSCI. 



Figure 6: Idem as the Figure 4. using the covariance matrix 
F5 instead of Fj . 



tested on simulated data and compared to the standard es- 
timation techniques that requires either 4 images or 2 inde- 
pendant images (OSCI). The results show that the proposed 
method gives good approximation of the degree of polariza- 
tion. This study needs to be extended with a theoritical anal- 
yses in order to precise the conditions of validity of the pro- 
posed approach. 
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Figure 5: Idem as the Figure 3. using the covariance matrix 
Fs instead of Fj . 



